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SUMMARY 

During the period of this report, January 15, 1969 to May 19 , 1971, 
work on this project was conducted in five major categories; In four of 
these categories, papers have been accepted for publication. 

In the first category, diffusion and permeation of C0 2 , 0 2 , N 2 , and 
N0 2 through poly tetra f luoroethy lene (PTFE ) were studied. Appendix A on 
this subject has been published in Macromolecules. This study indicates 
that the gases tested, except N0 2 , do not interact with the PTFE but dif- 
fuse through preexisting channels and voids. 

In a second category is the study of the diffusion, permeation, and 
solubility of simple gases (C0 2 , 0 2 , N 2 , CH 4 , C 2 H 6 , C 3 H 8 , and C 2 H 4 ) through 
a copolymer of hexaf luoropropylene and tetraf luoroethy lene (FEP) . The re- 
sults of this study are summarized in. Appendix B, which has been accepted 
for publication in Macromolecules. Existing theories of gas-polymer inter- 
actions are shown to be unable to explain satisfactorily the experimental 
data for FEP. 

In a third category, the viscous flow and diffusion of gases through 
small apertures were studied. A report on this study has been accepted 
for publication by the Journal of Applied Physics and the manuscript ap- 
pears as Appendix C. This study indicates that for the gases tested, C0 2 , . 
0 2 , and N 2 , the observed diffusion values are smaller than those calculated 
by. a nearly constant factor. 

A fourth category includes. a study of diffusion and permeation of 0 2 , 
N 2 , C0 2 , CH 4 , C 2 H 6 , and C 3 H 8 through nitroso rubber. Appendix D is the 
manuscript *of the paper, which has been accepted for publication by the 
Journal of Applied Polymer Science,. This study shows that the solubilities 
of these gases in nitroso rubber are 2 to 4 times those in natural rubber. 
The different behavior of C0 2 _ indicates a strong specific interaction. 
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Category 5 includes the results of gas transport studies on five 
materials submitted by JPL. This work is summarized in the following 
section . 

Appendix E lists other pubi „ions generated under this program or 
under a related NASA-JPL. program. 
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GAS TRANSPORT STUDIES ON MATERIALS SUBMITTED TO SRI BY JPL 


A . Carborane Siloxane 

The sample had a very rough, grainy surface; its thickness, which 
was measured with a micrometer, was nominally 22 mil. Some difficulties 
were encountered in sealing the membrane leak-tight into the cell but 
they were overcome by affixing on either side of the membrane 5-mil 
polyethylene gaskets with RTV silicone adhesive . 

The permeation and diffusion coefficients for C0 2 were determined 
at three temperatures. Hydrogen permeability was measured with nitrogen 
as carrier gas to achieve high detector sensitivity, and the detector 
was calibrated with a hydrogen-nitrogen mixture of known composition. 
After completion of the experimental series at maximum temperature 
(99. 5° C), the membrane was badly deformed , having assumed the shape of 
the cell, i.e., the material had flowed. 

The data for both COg and Hg are shown in Table 1. The permeation 
rate for hydrogen was so fast at higher temperatures that the calculated 
diffusion coefficients are not reliable; thus, the heat of diffusion 
quoted represents only a lower limit . 













B. 


Nitroso Rubber (CTVR-CTA cured batch 1254/1256) 


A nitroso rubber sample 1 mm thick was studied . Table 2 lists the 
permeation and diffusion coefficients of C0 2 measured at five tempera- 
tures . 


Table 2 

PERMEATION AND DIFFUSION COEFFICIENTS OF C0 2 
THROUGH A 1-mir. NITROSO RUBBER MEMBRANE 


t(°C) 

P x 10® 

/ cc(stp)cm ) 

D x 10 7 
(cm 2 /sec) 

\crci* sec cm.Hg/ 

25 

15.2 

6 .8 

30 

17.4 

8.6 

35 

19.1 

9.5 

40 

23.9 

11.4 

45 

26.1 

13.7 


The heats of permeation and diffusion calculated from these data are 
5.2 and 7.3 kcal/mole , respectively. 

Some data on permeation . through nitroso rubber are contained in a 
report by Thiokol Chemical Corp. (RMD 5088-5098-S , Contract No. NAS 
7-451). By making some guesses concerning the units used in that re- 
port, we calculate their permeation coefficient for CO^ at 25°C to agree 
with our value within 3%. 

In Table 3 the CC^ permeation data. for the nitroso rubber are com- 
pared with those of j>olyisoprene ana . of low density polyethylene; they 
are taken from the Polymer Handbook. 
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It is Seen that the nitroso rubber is very permeable, comparable 
to an amorphous hydrocarbon rubber. 

The interaction of N0 2 -N 2 0 4 with the rubber was explored also. Our 
source of the permeant was an N 2 0 4 cylinder which was kept at about 
30°C. Since the boiling point of at 1 a<tm is 21°C, a sufficient 

overpressure is generated to drive the vapor through the permeation 
cell. (The emerging N 2 0 4 . is either condensed in a dry ice trap or ab- 
sorbed in water.) The permeant passing into the He stream is present 
as pure N0 2 because of its low partial pressure (always well below 0.1%). 
The detector sensitivity for NO^ was assumed to be the same as that for 
CQ2 . 

One permeation run was carried out at 48.2°C.. From the steady state 

signal, which was reached after about 1-1/2 hr, a permeation coefficient 
_ 7 cc (stp)cm 

of 1.9 x 10 — 5 was calculated. A nominal diffusion coeffi- 

cm sec cm Hg 

cient of 7 x 10" 7 cm 2 /sec was derived, by our standard method. However, 
this number has at best only a semiquant itative significance, since two 
species, NOg and N 2 0 4 , are present. Thus the diffusion process is com- 
plex, and our simple procedure for determining D is questionable. 

When the sample was degassed , extensive tailing was. observed ; after 
the signal had dropped to about half its value, several hours were re- 
quired to reach th<* base line again. Thus sorption of N0 2 -N 2 0 4 by ni- 
troso rubber does not appear to.be reversible. The observed effects 
were shown to be due to the polymer itself, and not to slow sorption 
processes of N0 ? on the system walls: 1 to 50 (j,l air and , in paralr- 

lel experiments, N0 2 -Ns0 4 were injected with a microsyringe into the 
carrier gas stream. For both gases, the appearance, time of the signal 
was the same, the peaks Were symmetrical , and the peak heights were ap- 
proximately proportional to the amount injected. 
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c. 


Silicone Membrane (SR 6750-1) 


The hydrogen permeation was investigated over a temperature range 

of 50 .0 101° C . The data are shown in Table 4. The heats of permeation 

E and diffusion E were found to be fi =3.1 and E = 2.5 kcal/mole , 

P D P D 

and the permeation and diffusion coefficients at 50° C Were 

cc(stp)cm 

P = 5.2 X 10 8 — s — 

50 cm sec cm Hg 


and 


50 


4.5 x 10" 5 cm 2 /sec 


Table 4 

HYDROGEN PERMEATION DATA 
FOR SILICONE MATERIAL (SR 6750-1) 


t 

(°C) 

D X 10 s 
( cm 2 /sec) 

P X 10 8 

( cc(stp)cm \ 

S X 10 3 
/ cc(stp) \ 

\cm 3 cm Hg/ 
(calculated ) 

\cm 2 sec cm Hg/ 

101.1 

7.50 

9.94 

1.32 

91.7 

7.12 

8.90 

1.25 

80.1 

6.38 

7.74 

1.21. 

70.3 

5.68 

6.82 

1.20 

63.5 

5.37 

6.24 

1.16 

50,1 

4.54 

5.18 . 

1.14 


D . Krytox Coating on Teflon 

We explored whether a Krytox coating would reduce permeation rates 
through Teflon. A 4. 4-mil Teflon membrane, of the same origin as those 
used in the principal studies, was mounted in our permeation unit, and 
the permeation rate of C0^ was determined at five temperatures. The 
values agreed within 5% with the previously reported data. A 3-wt % 
solution of Krytox in Freon TF was then applied to the membrane in situ 
and the solvent evaporated. The average thickness of the residue was 
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about 0.5 mil. The coating whs not uniform because the membrane was 
not perfectly flat; visual 1 ispection indicated, however, that some ma- 
terial was deposited over the entire surface. The membrane was then 
carefully degassed and the permeation of CO^ was remeasured at two tem- 
peratures. The permeation rates were, within the precision of the meas- 
urements, identical with those for the coated membrane, 

E. Permeation and Diffusion of C0 2 and 0 2 Through FEP Coated Glass 
Cloth 

Permeation and diffusion measurements were carried out on samples of 
glass cloth inpregnated with FEP. The glass fibers had a diameter of 0.7 
mil and the open areas between the fibers were about one-half of the total 
area. The thickness of the glass-supported FEP membrane was calculated 
from its weight to be 5.5 mil. The directly ^measured thickness, which 
varied by about 20 $ since the glass cloth was not uniformly embedded in 
the plastic, was about 8 mil. The value 5.5 mil was used for the calcu- 
lations. 

Diffusion and permeation of C0 2 and 0 2 were measured at 65.0, 50.0, 
and 34 . 5°C . Reasonably straight Arrhenius plots were obtained. In Table 5, 
diffusion and permeation coefficients at 25°C, obtained by extrapolation 
of the Arrhenius plots, and the heat of diffusion and permeation of glass- 
supported FEP are compared with those of unsupported FEP. The agreement 
is surprisingly good if one considers the uncertainties in the present 
measurements. Thus the glass cloth may impart additional strength to the 
FEP membrane without changing its transport properties. 
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Table 5 

COMPARISON OF TRANSPORT PROPERTIES OF 
GLASS-CLOTH-SUPPORTED FEP(l) AND UNSUPPORTED FEP(2) 



C0 2 0.94 1.05 8.40 8.75 1.00 1.27 5.3 5.2 

0 2 1.50 1.84 8.40. 8.29 0.50 0.49 5.2 6.1 
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DIFFUSION AND PERMEATION OF 
OXYGEN, NITROGEN, CARBON DIOXIDE, AND 
NITROGEN DIOXIDE THROUGH POLYTETRAFLUOROETHYLENE 


by 

R. A. Pasternak, M. V. Christensen, and J. Heller 
Polymer Chemistry Program 
Stanford Research Institute 
Menlo Park, California 94025 


ABSTRACT 

The permeation of C0 2 , 0 2 , N 2 , and N 2 0 4 -N0 2 equilibrium mixtures 
through polytetraf luoroethy lene (PTFE) has b^. . investigated systemat- 
ically. PTFE membranes between 1.5 and 5.7 mil thick were studied over 
a temperature range of 30 to 116°C. For all the gases, the permeation 
and diffusion coefficients were found to be independent of membrane thick- 
ness, within the precision of the measurements. Linear Arrhenius plots 
were obtained for C0 2 , 0 2 , and N 2 . For the N 2 0 4 -N0 2 system, the temper- 
ature dependence of both permeation and diffusion was -complex; the per- 
meation data could be. explained quantitatively in terms of the N 2 0 4 -N0 2 
equilibrium. . 

The parameters describing the diffusion and solution process indicate 
that the C0 2 , 0 2 , N 2 , and N 2 0 4 molecules do not significantly interact 
with PTFE but move through preexisting channels and voids, whereas N0 2 
interacts strongly, though reversibly', with-the polymer. 


* Published in MacromdleculeS 3, 366 (1970). 
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DIFFUSION AND PERMEATION OF 
OXYGEN, NITROGEN, CARBON DIOXIDE, AND 
NITROGEN DIOXIDE THROUGH POLYTKTRAFLUOROETHYLENE 


R. A, Pasternak, M. V. Christensen, and J. Heller 
Polymer Chemistry Department 
Stanford Research Institute 
Menlo Park, California 94025 


Introduction 

The development of advanced liquid propulsion systems requires the 
use of expulsion bladders for storable propellants in the space environ- 
ment . The bladder material . must not only be chemically inert toward 
strong oxidants such as nitrogen dioxide but must also have a very low 
permeability toward these agents. Because of its chemical resistance and 
good mechanical properties, polytetrafluoroethylene (PTFE) has received 
serious consideration as a bladder material and several space vehicles 
have been launched using PTFE expulsion bladders . 

However, permeation of oxidants through PTFE, while not detrimental 
in short missions, is of serious, concern in long missions. Since only 
very limited data for the diffusion and permeation of noncondensible gases 
and none for nitrogen dioxide through PTFE are found in the literature, 1 * 2 
a systematic study of these gases was undertaken in an attempt to under- 
stand the details of the interaction between permeants ana PTFE. 

Experiment al 

The dynamic method for measuring permeation and diffusion rates of 
gases and vapors in polymer membranes has been described in detail else- 
where 3 and is here only summarized . The two compartments of the permea- 
tion cell, which are separated by the test membrane, are -open to the 
atmosphere. The permeant, gas or vapor, is admitted to one compartment, 
replacing a carrier gas, usually helium. A second carrier g as stre_sm 
flows at constant rate through the other compartment and sweeps the 


permeant, which diffuses through the membrane, to the thermal conductivity 
detector. The detector signal, which is originally zeroed with the car- 
rier gas in both compartments, is at any instant proportional to the. per- 
meation rate. When the permeant is again replaced by the carrier gas, 
the membrane is degassed, and the signal returns to the base line. Thus, 
repeat measurements can be made with the same membrane over a wide tem- 
perature range and for diverse permeants . 

The instrument used in this study is a prototype of the Polymer Per- 
meation Analyzer (PPA) .* The cell and the plumbing were stainless steel; 
the effective membrane area was 3 .9 cm 2 . The detector was a Loeneo^ ther- 
mal conductivity cell which was run at 140 mA and thermostated at 80°C . 

Its sensitivity for C0 2 , 144 ppm/mV, did not change significantly with time, 
even after exposure to N0 2 . The combination of the detector with a Sargent 
recorder (full scale deflection of 250 mm at an input of 0.4 mV) resulted 
in a detection limit of 0.2 ppm of C0 2 in helium. Helium was used as 
carrier gas for the study of all permeants, except of hydrogen, for which 
nitrogen was used . (With N 2 as carrier gas the detector current was 
50 mA.) The flow rates of the carrier gas were set between 0.5 and 1.5 
ml/min, the higher values being, employed at high diffusion rates. At 
steadj'-state permeation, strict inverse proportionality between signal 
and flow rate was found as expected for a detector with linear response. 

The detector sensitivity, which is a function of the gas or vapor 
sensed, was determined for C0 2 and H 2 0 . Helium with 0.21 vol $> C0 2 and 
nitrogen with 0.20 vol ^ H 2 , were used as test gases. Moreover, the per- 
meation rates of C0 2 , 0 2 , and N 2 through a PTFE membrane were measured in 
a Barrer-type permeation .cell (Linde Model CS 135), and. compared, with the 
equivalent signals obtained in .the PPA. The permeation coefficients for 
C0 2 derived by the two independent approaches differed by 7$; the data 
reported here are based on the calibration of the PPA . The mole response 
factors of the detector for 0 2 and N 2 (relative to C0 2 ) were found to be 
1 .20 and 1 .22 ; the factor for N0 2 was assumed to be unity because its 


* Manufactured by Dohrmann Instrument Co. (Infotronics) , Mt . View. 
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molecular weight differs by only two units from that of C0 2 . No N 2 0 4 
reached the detector since in the helium stream the highest partial pres- 
sure of the nitrogen oxides, expressed as N0 2 , was less than 0,1$. Final- 
ly, the equilibrium solubility of N 2 0 4 -N0 2 in PTFE at a total pressure of 
one atmosphere was measured over a limited temperature range with a quartz 
spring balance . 

High purity C0 2 , 0 2 , and N 2 were used without further purification 
or drying; the N 2 0 4 -N0 2 mixture was bled into the cell from a tank held 
at 30°C. The PTFE used in this study had a density of 2.18 g/cc at 18°C 
as determined by immersion in water. Membranes were cut with a microtome 
from a solid polymer block of Dupont Teflon, commercial stock, and their 
thicknesses were determined by a micrometer with a reading precision of 
0.05 mil. Only membranes with thickness variation of less than 5$ were 
used. 


Evaluation of Permeation and Diffusion Coefficients 

In Fig. la, a typical permeation curve for C0 2 is shown. The plateau 
represents steady-state permeation. The permeation coefficient P is cal- 
culated by the formula 


P 


f .q«k«a •Sqq.'L 
A *p 


("cm 5 ' • 


cc»cm 
sec*cm Hg, 


where f is the flow rate, cr the detector sensitivity, k the Molecular 
Response . Factor , a the detector attenuation, the steady-state signal, 
t the membrane thickness, A the membrane area, and p the upstream permeant 
pressure . The permeant pressure is taken as 76 cm Hg although it varied 
slightly. 

The steady-state data for N 2 0 4 -N0 2 are reported as the measured flux 
of N0 2 at a total pressure of one atmosphere, normalized to unit area and 
thickness. Formalistically , the flux is identical to a nominal permeation 
coefficient at atmospheric pressure . 

The diffusion coefficient was derived from the ascending "and descend- 
ing- branches of Fig. la which represent the transient permeation rates. 
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The solution of the differential equation describing the transient state 
has been discussed in detail elsewhere. 3 It is to a good approximation 



X *exp (-X 2 ) 


where X 2 = l 2 /4Dt , S is the signal at time t, and D is the concentration- 
independent diffusion coefficient. This equation describes both sorption 
and degassing, which are identical except for reversal in the direction 
of the signal. A plot of S/S^ vs l/x 2 is the theoretical, normalized, 
permeation rate curve, if, for the same values of S/S^, the normalized 
times l/X 2 are plotted vs the experimental times t, the straight line 
through the origin should have a slope equal to 4D/£ 2 . Figure lb repre- 
sents this evaluation of the ascending and descending branches of Fig. la 
The points lie on a straight line, and close agreement is found between 
the two branches. However, a positive displacement on the time axis is 
found, because a finite time is required for the permeant to reach the 
upper compartment, and for the permeated gas to reach the detector. 

A simplified method has been adopted for the evaluation of the large 

a 1 /x 2 

amount of data. The slope, which is given by the ratio — ^ — , was eval- 
uated from two sets of two experimental points each, namely for S/S w = 
0.1, 0.4 and S/S^ =0.4, 0.7. For the noncondensible gases, the two 
values agreed within 5 and were averaged . 

For the N0 2 -N 2 0 4 system the determination .of diffusion coefficients 
by the described formalism becomes questionable. The simultaneous diffu- 
sion of two interacting components must result in a complex rate law; 
unfortunately, the differential equation cannot be solved explicitly in 
this case * The experimental curves have shapes similar to those for the 
noncondensible gases (Fig. 2a); however, plots of l/x 2 vs t (Fig. 2b), 
though approximately linear, show a slight curvature and they intercept 
the time axes at negative values. Both effects diminish with increasing 
temperature, as would be expected from a consideration of the N 2 0 4 -N0 2 
equilibrium. The nominal diffusion coefficients reported here were some- 
what arbitrarily derived from the slopes for the interval S/S^ =0.1, 0.4 


Results and Discussion 


1 , Carbon Dioxide , Oxygen, and Nitrogen 

In Fig. 3 the logarithm of the permeation coefficients for C0 2 , 0 2 , 
and N 2 derived from the data are plotted versus l/T; Fig. 4 is an equiv- 
alent plot of the diffusion coefficients. For C0 2 , membranes of thickness 
1,5, 3.2, 3.5, 4.3, and 5.7 mil were used, but no diffusion coefficients 
are given for the thinnest membrane, since the transient permeation rate 
was too fast to yield reliable values. The 1.5-mil membrane was omitted 
in the 0 2 and N 2 studies. 

The complete data cover a temperature range of 30 to 116°C; however, 
the individual membranes were studied over smaller intervals, the range 
depending on the membrane thickness. Also, for any given membrane, per- 
meation coefficients are usually given to higher temperatures than diffu- 
sion coefficients because the determination of the latter becomes unreli- 
able for fast transient rates. 

Since the agreement between the data for different membranes was very 
good, they were all combined in a least square treatment. The straight 
lines in the figures represent the results. The mean square root devia- 
tion of the individual data points from the best line is between 3 and 5$ 
for all plots. Deviations for an individual membrane do not appear to be 
systematic or significant, and are well within the range of uncertainties 
originating from thickness measurement and slight variacion in flow rate 
of the carrier gas . 

In Table 1 the diffusion and permeation coefficients at 25°C, D 25 and 
P 25 , the preexponential factors D 0 and P 0 , and the. heats E^ and E p , are 
summarized. In addition, the solubility coefficients calculated from 
s = P/D are listed. The permeation coefficient and the heat of permea- 
tion of H 2 , which was measured for one membrane only and at four temper- 
atures between 24 and 48°C, are included. Permeation coefficients for 
the four gases in FEP given in the literature 2 are shown also. 

The diffusion and permeation coefficients of C0 2 , 0 2 , and N 2 , -plotted 
in Figs. 3 and 4, confirm convincingly that the generally assumed simple 




diffusion law holds; the independence of the coefficients oh membrane 
thickness, in particular, proves that the sorption process on the membrane 
surface is fast and is not the rate determining step, the data also show 
that the activation energies of diffusion and permeation are independent 
of temperature, over the rather wide temperature range of this study. 

The activation energies of diffusion E^ of the three gases (Table 1) 
are quite similar. However, in view of the high precision of the data, 
the differences, though small, are significant. This is also indicated 
by the trend in the preexponential factors; their logarithms, i .e . , the 
entropies of activation, are approximately a linear function of the acti- 
vation energies of diffusion, as has been found in other systems. 1 The 
trend in E^ appears to be associated with the molecular sizes of the gases . 
Oxygen has a significantly smaller Van der Waal's diameter (2 .92A) than 
N 2 (3 .15A) and C0 2 (3 .23A) , 4 Moreover, since the cross section perpendic- 
ular to the molecular axis is smaller for C0 2 than for N* , the lower E 
value for C0 2 may indicate that the molecules permeate by preference in 
the direction of their molecular axis . 

The heats of solution of the three gases in PTFE and their heats of 
condensation (last column) are compared in Table 2. The close numerical 
agreement is probably fortuitous; nevertheless it indicates that solution 
in PTFE is similar to a condensation process with minimal interaction be- 
tween the polymer and the solute . 

A comparison of diffusion and solubility data of PTFE with those of 
the structurally similar linear polyethylene is of some interest (Table 2). 
The polyethylene data are for a Phillips type 22$ amorphous polymer and 
are taken from two papers by Michaels and Bixler 5 ; they are converted to 
the units used here . The solubilities refer. to the amorphous phase and 
are calculated assuming the solubility in the crystalline phase to. be 
negligible. The PTFE used was 40$’ amorphous as determined from its den- 
sity . 

The activation energy of diffusion of the three gases in PTFE is about 
two kcal/mole lower than in polyethylene. This is surprising, since PTFE 
has very stiff backbone chains of high viscosity. Apparently, as proposed 


Table 2 


DIFFUSION AND SOLUBILITY DATA OF PTFE (T) 
AND OF LINEAR POLYETHYLENE (PE)* 


Gas 

E 

D 

S 25 

X 10 2 

. * H S 

AH , 

cond 

T 

PE 

T 

PE 

T 

PE 

co 2 

6 .84 

8 .5 

3 .1 

0.60 

- 3 . 5.0 

-1 .3 

-4.0 + 

0* 

6 .28 

8 .8 

0 .69 

0.10 

- 1 .73 

-0 .4 

-1 .6 

n 2 

7 .12 

9 .0 

0.39 

0.054 

-1 .30 

0.5 

-1 .3 


Note; AH in kcal/mole; s 2S in — 3 — 

cm cm Hg 

* Kef. 5. 

t At triple point, -56°C. 


originally by Brandt and Anysas, 1 diffusion of small molecules in. this 
polymer proceeds through preformed channels and cavities and requires lit- 
tle movement of the polymer chains. 

The higher gas solubilities at 25°C, more than six times those in 
polyethylene, and the higher heats of solution may again indicate poor 
packing of the PTFE chains in the amorphous regions. 


2 . Nitrogen Dioxide 

Figure. 5 summarizes the data for the N 2 0 4 -N0 2 system obtained with 
the same five membranes as for C0 2 . Repeat runs were reproducible r;o n 
after long time exposure of the membranes to the .strong oxidant. More- 
over, ..data for C0 2 obtained before and after exposure to N0 2 were identi- 
cal within the precision of the measurements, evidence that PTFE is not 
attacked by N0 2 . The nominal diffusion coefficients (see. preceding sec- 
tion) lie on a straight line (curve 1) except./at the lower temperatures. 
The steady-state permeation data closely fit the theoretical curve 2 (see. 
below),' except for some points for the 4.3-mil membrane at lower tempera- 
tures. 
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Finally, in Table 3 the solubilities of the N 2 0 4 -N0 2 equilibrium mix- 
tures at a total pressure of one atmosphere are given for two PTFE sam- 
ples ; the agreement is quite satisfactory. 


Table 3 


SOLUBILITY IN PTFE OF N 2 0 4 -N0 2 
EQUILIBRIUM MIXTURE AT 1 ATM PRESSURE 


T 

(°C) 

s(g/cc Atm) 

s(cc N0j (STP) /cc Atm) 
(average) 

Sample 1 

Sample 2 

30.1 

0 .039 


19.0 

35 .4 

0 .027 


13 .0 

40.6 

0.019 


9.5 


The close agreement of the apparent diffusion coefficients and the 
normalized permeation rates for different membrane thicknesses, shown in 
Fig. 5, indicate again that the adsorption process is fast, £nd that the 
transport within, the membrane is the rate determining step. This trans- 
port mechanism however is complex, not only because two species diffuse, 
but also because the relative concentrations of N 2 0 4 and N0 2 in the gas 
phase and in the polymer are functions of temperature and total pressure 
(or concentration) . 


The equilibrium in the gas phase is given by 6 


(P 

~P 


) 2 

NO* 

n 2 o 4 



7.1 


10 


9 -14,600 A 

e>XP RT~ Atm 


( 1 ) 


At a. total' pressure of one atmosphere on the upstream side of the membrane 


(our standard experimental condition) the partial pressure of NQ 2 is 0.36 
atm at 30 C, and 0.99 atm at 116°C . (On the downstream side, virtually 
pure N0 2 is released since the total pressure of permeant in the helium 
stream was in all experiments -below 10~ 3 atm.) 
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Since at the highest temperatures of this study the gas phase con- 
sists of virtually pure N0 2 , it is reasonable to assume that the measured 
nominal diffusion and permeation coefficients are the true coefficients 
for N0 2 . Actually, the diffusion coefficient curve is at the higher tem- 
peratures linear over a much wider interval than anticipated. The diffu- 
sion parameters derived from it are ascribed to N0 2 and listed in Table 4 


Table 4 

DIFFUSION, PERMEATION, AND SOLUBILITY DATA 
FOR N0 2 AND N 2 0 4 IN PTFE 
(For Units and Definitions, see Table 2) 


Gas 

D 25 xio 7 

Do 

e d 

P 2S X 10 9 

P o X10 6 


s 2S X10 2 

s o X10 9 

AH 

s 

no 2 

0 .037 

70 

14 .0 

1 .57 

0 .116 

2 .55 

43 

1 .66 

-11 .4 

^2^4 

0.25 

— 

— 

3 .75 

0 .0013 

-0 .5 

15 


— 


In our analysis of the observed nonlinear temperature dependence of 
permeation, we make the following postulates: (1) the solubility of either 

component obeys Henry’s Law 


c ;_ = s ■ h < 2 > 

where C° is the saturation concentration in the polymer, s_ is the solu- 
bility coefficient, and p i is the pressure in the gas phase; (2) both the 
diffusion and solubility coefficients are independent of concentration. 
This postulate probably holds only approximately since the total equilib- 
rium amount of N0 2 + N 2 0 4 dissolved at the lowest temperature of the ex- 
perimental series is about 2$ by weight ; (3) the equilibrium N 2 0 4 ^ 2N0 2 
is maintained at all times in the gas phase and everywhere within the mem- 
brane, even during transient permeation. 
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It follows from 
where in the polymer 


the three postulates that the concentrations any- 
satisfy, at all times, the equilibrium condition 



Diffusion through the membrane is described by the two differential 
equations 


dC 


NO, 


d 2 C. 


dt 


-D. 


NO, 


NO, 


dC 


N ^l = _ D . 


d 2 C. 


N,0. 


dt 


N 2 0 4 


dx 2 


* - H 


(4a) 


(4b) 


R is the net rate of conversion of N 2 0 4 to N0 2 ; the factor 2 in Eq . (4a) 
accounts for the production of two N0 2 molecules for every dissociating 
N 2 0 4 molecule . 

Equations (4a) and (4b) are not independent, but are related by Eq . 

(3) . By combining these three equations, a differential equation in , 

NU 2 

time t, and position x can be obtained* Unfortunately, this equation can- 
not be solved except numerically, and therefore the diffusion coefficients 
cannot be derived explicitly from the experimental curves. However, a 
solution can be obtained for the steady state, when, the concentrations of 
both N0 2 and N 2 0 4 anywhere in the membrane are independent of time. The 
flux F of permeant through the membrane of thickness l 


dC 


F = -D. 


NO, 


dC. 


NO, dx 


- 2D, 


N,0, 


n 2 o 4 


dx 


2^4 


(5) 


is then constant, i .e . , also independent of x. The factor 2 has been 
introduced to express the flux in terms of the N0 2 species, which appears 
exclusively on the downstream side of the membrane. 
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At steady state the boundary conditions are 


at x = 0 


C - r° 
N0 2 "NO 2 


C — c° 
n 2 o 4 n 2 o 4 


at x = t 


C = C =0 

no 2 n 2 o 4 


Integration of Eq. (5) between these limits yields 


F * / = D C° + 2D C° 

^ N0 2 N0 2 n 2 o 4 n 2 o 4 


( 6 ) 


On introducing Henry's Law, Eq . (1), we obtain 


F ' 1 d no 2 s no 2 p no 2 + 2D n 2 o 4 s n 2 o 4 p n 2 o 4 


(7) 


We define: the flux at unit thickness F 0 = F*£, and 


the permeation coefficients P = D s : 

N0 2 N0 2 N0 2 ’ 


p n 2 o 4 d n 2 o 4 s n 2 o 4 


Since in our experiments the pressure condition holds 


p no 2 + p n 2 o 4 = 1 atm ' 


the flux equation assumes the simple form 


F ° _ p no 2 p no 2 J" 2(1_p no 2 )p n 2 o 4 


( 8 ) 
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Formalistically this equation is identical with the expression that de- 
scribes the permeation of two independent components; however, in the 
N0 2 -N 2 0 4 case, the concentrations of the two species in the membrane are 
linear functions of x but are related by the' mass law. 

It can be easily seen that the permeation data plotted in Fig. 5 are 
in the form of this normalized flux F 0 ; they can be analyzed in terms of 
Eq. (8). We postulate that both P^ Q and P N Q show an exponential depend 
ence on the reciprocal temperature. Since at the highest temperatures of 
the experiments the gas phase is virtually pure N0 2 , the permeation coef- 
ficients P N0 ^ are given by the asymptote to the high temperature data, 
curve 3 in Fig. 5, and can be extrapolated to lower temperatures. More- 
over, the partial pressures of p NQ are given by equilibrium Eq . (1). 
Consequently, only the permeation coefficient P N of Eq . (8) is unknown 
and can be calculated at different temperatures from the smoothed, experi- 
mental F 0 values. After minor adjustments in the slope of the P curve 
, . r N0 2 ’ 

a linear P^ curve was obtained [(4) in Fig. 5).}. Then, as a check, the 

total flux F 0 was recalculated from the two permeation curves (3) and (4) . 
The agreement of the calculated curve (2) with the experimental data is 
excellent , except for some points derived from the 4 .3-mil membrane . 

The parameters derived from the two permeation curves (3) and (4) are 
listed in Table 4. The. coefficients, which in Fig. 5 are given in atmos- 
-pheric pressure units,. are converted to cm Hg units to facilitate compari- 
son with the C0 2 , 0 2 , and N 2 data. Approximate solubility and diffusion 
coefficients of N 2 0 4 at 25°C are given also. The solubility coefficient 
was derived from the directly measured solubilities (Table 2) by extrapo- 
lation to 25°C. A correction was made for the solubility of the N0 2 pres- 
ent in the gas phase at 25°C . 

A comparison of the solubility and diffusion coefficients of N 2 0 4 and 
of C 0 2 , 0 2 , N 2 indicates a common- pattern . The Solubility coefficients 
of the gases correlate quite well with their boiling points, as has been 
found for other gas-polymer systems. 7 Similarly, the reduction in the 
N 2 0 4 diffusion coefficient, by a factor of about 4 compared to C0 2 appears 
to be reasonable in view of the much larger size of the 1^0, molecule. 

Thus, N 2 0 4 , -like the other gases, apparently does not interact With PTFE. 
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In contrast, :he N0 2 data do not fit the pattern at all. Although 
the N0 2 molecule has approximately the same size as C0 2 , its activation 
energy of diffusion is about twice and its heat of solution about three 
times that of C0 2 . Also, its diffusion and solubility coefficients at 25°C 
are smaller and larger, respectively by orders 6 f magnitude. These obser- 
vations suggest that the highly polar N0 2 molecule interacts strongly with 
PTFE . 

A systematic analysis of the data in terms of molecular dimensions 
and interchain forces is not only of considerable interest, but is essen- 
tial for a more detailed and quantitative understanding of diffusion mech- 
anisms . This analysis must however be postponed until data for a larger 
number of permeants and PTFE samples of different degree of crystallinity 
are available. 
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ILLUSTRATIONS 


1 (a) Experimental curve for permeation of C0 2 through a 4 ,3-miI 
PTFE membrane at 90.6°C. Arbitrary units. C0 2 was intro- 
duced at A; He at B. 

(b) Plot of l/x 2 vs time t at S/S^ = 0 .1 f 0 .2 , ... ,0 .9 . 

^Ascending branch. □ Descending branch. 

2. (a) Experimental curve for permeation of N 2 0 4 -N0 2 through a 
4.3-mil PTFE membrane at 49.7°C. Arbitrary units. 

(b) Plot of l/x 2 vs experimental time t at s/s = 0.1, 0.2,..,, 0.8. 

3 Arrhenius plots of permeation coefficients in PTFE. 

Arrhenius plots of diffusion coefficients in PTFE. The C0 2 
curve is shifted by a factor of 2 . 

Permeation and diffusion coefficients of N 2 0 4 -N0 2 equilibrium 
mixtures in PTFE . 

(1) Nominal permeation coefficients of mixture . 

(2) Nominal diffusion coefficients of mixture. 

(3) Permeation coefficients of N0 2 . 

(4) Permeation coefficients of N 2 0 4 . 
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DIFFUSION AND SOLUBILITY OF SIMPLE GASES THROUGH A 
COPOLYMER OF HEXAFLUQRQPRQPYLENE AND TETRAFLUOROETHYLENE 

by 

R. A. Pasternak, G. L. Burns, and J. Heller 
Polymer Chemistry Program 
Stanford Research Institute 
Menlo Park, California 94025 


ABSTRACT 


The transport of seven gases , N 2 , 0 2 , C0 2 , CK 4 , C 2 H 6 , C 3 1I 8 , and C 2 H 4 
through FEP , a copolymer of hexaf luoropropy lene and tetraf luoroethy lene , 
has been studied systematically over a temperature range of about 25 to 
85°C. For all the gases studied, linear Arrhenius plots of both permea- 
tion and diffusion coefficients were obtained. Empirical linear correla- 
tions both between the logarithm of the solubility at 25° and the heat of 
solution, and the Lennard-Jones force constant of the gases were found. 
Furthermore, the standard entropies of solution and of diffusion were lin 
early dependent on the associated enthalpies. A comparison of the FEP 
data with those for polyethylene taken from the literature points up the 
limitations of existing theories of transport through polymers.. Finally, 
the comparison with data on poly tetraf luoroethy lene (PTFE) indicates that 
the properties of the latter may be due at least in part, to micropores 
or grain boundaries. 


* To be published ini Macromoiecules . 
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DIFFUSION AND SOLUBILITY OF SIMPLE GASES THROUGH 
A COPOLYMER OF HEXAFLUOROPROPYLENE AND TETRAFLUOROETHYLENE 


By R . A. Pasternak, G. L. Burns, and J. Heller 
Stanford Research Institute 
Menlo Park, California 94025 


In a recent study of gas transport in poly tetraf luoroethylene (PTFE ) 
it was found that the activation energies of diffusion of N 2 , 0 2 , and C0 2 
were approximately 2 kcal/mole lower than they were in polyethylene (PE). 1 

.J£ince_complete fluorination of polyolefins leads to considerable stiffen- 
ing of the carbon backbone , this is an unexpected result which suggests 
that diffusion progresses partly through micropores in the PTFE. Such an 
explanation appears reasonable since PTFE of sufficiently high molecular 
weight to be usable does not flow in the melt above its crystalline melt- 
ing point of 327°C ; fabrication requires compacting PTFE powder under high 
pressure near room temperature and then sintering at temperatures around 
380°C. The final product usually obtained by machining very likely con- 
tains micropores. 

The present study of diffusion, through a copolymer of hexaf luoropro- 
pylene and tetraf luoroethylene (FEP) was undertaken partly to verify this 
postulated diffusion mechanism. FEP has a crystalline melting point near 
290°C, and whereas it. retains most of the properties of PTFE\ it has a 
melt viscosity in a range which allows fabrication by conventional tech- 
niques. Thus, FEP films can be expected . to be free of micropores. 

Experimental Section 

Measuring Method 

The dynamicL 4 nethod for measuring permeation and diffusion rates of 
gases and vapors_ in polymery membranes and the evaluation of the perinea- . 
tion and diffusion coefficients from the raw data have been described 
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in detail elsewhere . 2 » 3 The instrument used here was an engineering model 
of the Polymer Permeation Analyzer.* This instrument is equipped with 
three permeation cells, which can accept membranes 0 f different thickness, 
and thus it permits a wide choice of measuring conditions. A low volume 
valve permits selection of any of the three cells for study. 

The instrument is well suited to rapid study of a series of gases 
over a wide temperature range. In actual operation, a gas at atmospheric 
pressure is admitted from a manifold to the upstream side of the thermo- 
stated cells; after completion of the measurements, the manifold is evacu- 
ated with a mechanical pump, and another gas is admitted and the procedure 
repeated. 

As a test of reproducibility, the permeation rates of C0 2 at 35° 
through one piece of FEP mounted across all three cells were measured. 

They were found to agree within 2<f 0 , which is equal to the reproducibility 
of repeat measurements for one cell. 

The concentration of permeant gas in the helium carrier gas sweeping 
past the membrane was measured with a katharometer .. The absolute detector 
sensitivity for C0 2 was determined by means of a C0 2 -He mixture of known 
concentration. The sensitivities for the other gases studied relative to 
C0 2 were determined directly by a method which will be described elsewhere. 4 

The measured sensitivities, which compare quite well with those given 
in th> literature, 5 are: 

C0 2> 1; 0 2 , 1.21; N 2 , 1.14; CH 4 , 1.32; C 2 H 6 , 0.97; C 2 H 4 , 0.98; 

C 3 H 8 , 0.80. 

G ases . . — - 

High purity bottled gases were used without drying. Ethane and ethyl- 
ene were better than 99,0^, the other gases were 99.5 to 99.9^ pure. It 
should be recognized that contamination of a slowly - permeating gas with a 

* 

Manufactured by Infotronics Instrument Corn, , Mountain View, California. 
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rapidly permeating gas can greatly influence the data* For example, the 
permeation coefficient of C 3 Ii 8 at 25° is about 1$ of that of C0 2 . Thus, 
a Ifa contamination with C0 2 would nominally double the steady-state 
signal of C 3 H 8 , However, since a parallelism exists between permeation 
and diffusion coefficients, C0 2 appears on the downstream side long 
before C 3 H 8 , and two well-defined steps would be observed in the dynamic 
method used here. Indeed, a very low step preceding the main step was 
occassionally observed and was ascribed to contamination from the manifold. 

Membranes 

FEP membranes of different thickness obtained directly from DuPont 
and from local suppliers were compared at 35°C. Variations were found 
in the permeation coefficient for C0 2 of more than 1CK, larger than the 
uncertainty in thickness measurements (0.05 mil). Three membranes were 
finally selected which had average thicknesses of 1.10, 1,87, _and 5.26 
mils. The permeation coefficients for the first two membranes were identi- 
cal within the precision of the measurements, but that for the thick mem- 
brane was lower by about 10$. This discrepancy was found at all tempera- 
tures and for all gases, whenever a comparison could be made; the diffusion 
coefficients, however, did not appear to vary for the three membranes. 

The permeation data reported here are a 1.1 derived from the l.lQ-mil mem- 
brane; the diffusion data arc combined from all membranes since, depending 
on the gas and the temperature, steady-state permeation might be reached 
too fast for the thin membranes, and too slow for the thick membranes to 
determine reliably the diffusion coefficients. In the present instrument, 
the time to reach half the steady-state signal must be at least 30 sec to 
allow evaluation of the curves in terms of diffusion coefficients. 

The membranes were annealed for 24 hr at 95°C. No significant change 
in thickness was observed; however, both the permeation .and diffusion co- 
efficients at 35° had increased by _ about 10$. Occasional temperature 
cycling, between higher temperatures and 35° in the course of the study 
did not Show further changes in the transport parameters. 


Crystallinity 


Since no data on the crystallinity of FEP were given by the manufac- 
turer or in the literature f x-ray diffractometer spectra were taken for 
both FEP and PTFE . A few distinct lines were found superposed on diffuse 
bands which were identical for the two polymers; they are given in Table 1. 
Other lines are indicated but are barely distinguishable from the back- 
ground . 


Table 1 

DIFFRACTION DATA FOR PTFE AND FEP, 
INDEXED ACCORDING TO A HEXAGONAL LATTICE 



PTFE 

FEP 

hkl 

Peak Heights 

d 

a = b 

Peak Heights 

d 

a = b 


(arb units) 

(A) 

cb 

(arb units) 

d) 

(A) 

100 

356 

4.90 

5 . 65 

208 

5.01 

5.77 

110 

16 

2.83 

5.66 

5 

2.885 

5.77 

200 

( 20 

2.43 

5.58 

6 

2.49 

5,74 


( 20 

2.18 

— 

— 

— 

— 


Note: hkl are the indices, d the spacings, and a = b the 

hexagonal axes. 

The first three lines listed fit well a hexagonal lattice, with axes 
of 5.63 and 5.76 A for PTFE and FEP, respectively. The former value is 

o 

close to 5.66 A, reported previously . 6 ’ 7 

The close agreement in unit cell dimensions and the similar relative 
intensities of the lines suggest closely similar packing of the chains in 
the two polymers; apparently the bulky CF 3 groups of FEP can be accommo- 
dated in the crystal structure without increasing significantly the unit 
cell dimensions. 

The lines of FEP were found to be significantly wider than those of 
PTFE. The crystallite sizes (in the a-b plane) estimated, from the line 
Width were about 200 A. for FEP and at least twice as' large- for PTFE. 
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to summarize, the FEP contains a significant fraction of a highly 
dispersed crystalline phase, though probably less than PTFE. However, 
the present data are not sufficient to make a quantitative estimate of 
this fraction. 

Result S 

In Figs. 1 and 2 the logarithms of the permeation coefficients and 
of the diffusion coefficients, respectively, are plotted versus the re- 
ciprocal absolute temperatures. A temperature range of about 60° is cov- 
ered by the data; except the steady-state permeation rates for C 3 H 8 at 
the lowest two temperatures were too small to permit precise measurements. 
The straight lines drawn were obtained by least squares treatment of the 
data. The standard deviation of the data points from this line is about 
2 <joi and the uncertainty in the slopes is of about the same magnitude; no 
systematic deviations are observed. However, both permeation and diffu- 
sion coefficients measured, at 95° C (not shown in the figures) lie signi- 
ficantly above the straight lines; these deviations were shown to be re- 
versible by temperature cycling. We suspect a phase transition, which 
is also indicated by a discontinuity in the tensile modulus at about 80°C. 8 

In Table 2 the permeation and diffusion coefficients P and D at 25° 
and the associated preexponential factors and energies are listed, they 
are derived from the straight line plots according to the equations 


P - P Q exp(-AEp/RT) 
D = D q exp(-AE D /RT) 


( 1 ) 


The solubilities and parameters associated with them, which are also 
shown, are calculated from the relationship 


k = f = k Q exp(- AHg/RT) 


( 2 ) 


The uncertainties in AE p and AE p are estimated to be about 0.2 kcal, 
and in AH q about 0.3 kcal. Uncertainties in the preexponential factors 
are quite large because their determination involves a linear extrapolation 
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over a wide interval in l/T. For an avei'age temperature of 55°C (T = 

323°), an uncertainty of 0.1 kcal in the heats results in an uncertainty 
of 15$ in tie preexponential factor. 

Discussion 

Permeation 

In Table 3, the experimental permeation coefficients P 25 at 25°C and 
the heats of permeation E p of FEP (and PTFE ) are compared with those of a 
linear and of a branched polyethylene PE Grex (Grace) and Alathon 14 (Du- 
Pont) , 22$ and 59$ amorphous, respectively. The PE data are taken from 
Michaels and Bixler. 9 N 2 and 0 2 permeate faster through the polyfluoro- 
carbons than even through the low density PE. However, the relative rates 
shift in favor of PE with increasing boiling temperature (or with increas- 
ing molecular size) of the reactant. 

In the following discussion of the two basic processes, solubility 
and diffusion, which together control permeation, only the data for linear 
PE (Grex) are used for comparison, since the degree of crystallinity of 
this PE does not change with temperature. We postulate that the crystal- 
linity of the high-melting FEP, after annealing at 95°C, is constant also. 

Solubility 

The solubilities of gases in some nonpolar solvents at 25°C and atmos- 
pheric pressure have been found to be an exponential function of the Lennard- 
Jones force constant e/k , 10 or the boiling temperature T^. 1 1 The same func- 
tional relationship has been observed for the .solubilities of gases in PE, 9 
natural rubber, 12 and polyethylene terephthalate. 1 3 The heat of solution, 
in PE is a linear function, of e/k also. 9 Since both e/k and T^ are measures 
for the van der Waals 1 interaction forces of gases and differ only by an 
approximately constant factor, the discussion will be Limited to e/k. 

In Fig. 3 the logarithms 6f the observed gas Solubilities k 25 at 25°C, 
and in Fig. 4*the heats of solution AHg , are plotted versus e/k for FEP, 

The values of e/k are those given by Michaels and Bixler, except that 
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r- 





for C 2 H 4 , which is taken from Hildebrand et al . 10 The fit to straight 
lines is excellent (except for C0 2 which possibly interacts with FEP) , 
but may be somewhat fortuitous since the values of the force constant 
s/k are not quite certain . 11 The equations of the straight lines are 

log k 2S = 1.38 + 0.0052 e/k ( 3 ) 

AH g = -0.85 - 0.012 e/k (kcal mole" 1 ) ( 4 ) 


For PE the equivalent log k 25 and AH^ plots have slopes of 0.0095 and 
0.0164. 9 

The solubility k and the heat of solution AH are related by the 
expression 


log k = log k - 
o 


2.3 RT 


If k Q (or the standard entropy of solution, AS q = 2.3 log k /R) of the 
different solutions were the same, the log k versus AH g plot would be lin- 
ear with a slope of -i/2.3 RT = -0.73 kcal " 1 mole at 298°K. The solubil- 
ities in PE satisfy this correlation . 9 For FEP (see Table 2) the k is, 
however, not constant; it is within the precision of measurements a linear 
function of £H g (Fig 5) with a slope of -0.29 kcal-mole" 1 . 

The linear interdependence of the standard entropy and enthalpy of 
solution appears at least qualitatively plausible, since the higher heats 
of solution, i.e., stronger interactions, should be associated with re- 
duced degree of freedom of the dissolved molecules. 

A theoretical interpretation of the observed empirical correlations 
involving the Lennard-Jones force constant has been attempted by Michaels 
and Bixler . 9 It was based on the Hildebrand-Scatchard 1 4 model for the 
mixing of two liquids and that model’s adaptation for polymers by Flory- 
Huggins . 15 In addition, the-chemical potential of the gases was assumed 
to obey the simple Clausius-Clapeyron equation. A linear correlation 
between AH g and e/k- was derived with a slope closely matching the exper- 
imental one, (0.0156 versus 0.0164 for PE ). 9 This agreement between 
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observation and model is probably fortuitous because of the rather tenuous 
assumptions made for the model: (1) the TrOuton rule holds, i.e., the 

entropy of vaporization at the boiling point is equal to 0.020 keal/deg 
for all gases studied, (2) the (fictitious) heat of condensation is cOn- 
stant even above the critical point; and (3) the heat of mixing is always 
small relative to the heat of condensation, or is constant. Postulate 
(3), made implicitly, conflicts with the theory of mixing of liquids as 
formulated by Hildebrand and others 14 and with experimental data for such 
systems. The contradictions with theory are illustrated by the data sum- 
marized in Table 4. Indeed, for PE the heats of solution are by an approxi- 
mately constant amount more positive than the heats of condensation; in 
contrast, for FEP the difference is negative for the low-boiling gases but 
decreases with increasing boiling point and becomes positive for propane. 
However, according to the simple theory of solution the heat of mixing 
can only be positive. Even larger negative heats of solution were found 
for gases in glassy polyethylene terephthalate; 13 they were ascribed to 
strong exothermic adsorption of the gases on the walls of internal voids. 

The gas solubilities at 25° C are listed also in Table 4; for PE and 
PTFE they are corrected for crystallinity. The uncorrected solubilities 
in FEP which represent a lower limit are significantly higher than in PE 
for the low boiling gases; but the reverse holds for the higher boiling 
gases. 

The solubilities of N 2 , Q 2> and C0 2 „ in PTFE appear to be significantly 
higher than in FEP, but- the heats of solution are somewhat less negative. 
Both these differences may be explained by quasi-condensation of the gases 
in micropores of PTFE. 

Diffusion 

It is seen from Table 2 that the preexponential factor D q (.which is 
about proportional to the entropy of. -activation) varies significantly * 
between solutes. In Fig, 6 the logarithm of D q is plotted Versus the 
activation energy of diffusion E^ of the gases in FEP. The data can be 


table 4 


SOLUBILITY DATA FOR DIVERSE GASES IN 
PE f FEP , AND PTFE 


Ga s 

e/k 

5 

AH. 

b 

PE 2 

FEP 

PTFE 2 

Cbnd 4 

PE 4 

FEP 

PTFE 3 

n 2 

95 

0.041 

0.13 

0.30 

-1.33 

■ / " 

-1.89 

-1.30 

0 2 

118 

0.077 

0.21 

0.52 

-1.63 


-2.20 

-1.73 

C0 2 

189 

0.451 

0.92 

2.35 

-3.46 

^^9 

-3.75 

-3.50 

ch 4 

148 

0.203 

0.22 

— 

-1.96 

-0.7 

-2.58 

— 

C 2 H s 

243 

1.28 

0.71 

— 

-3.51 

-2.3 

-3.61 

— 

c 3 h 8 

284 

3.97 

1.40 

— 

-4.49 

-3.0 

-4.08 

— 

c 2 k 4 

200 


0.49 


*3.24 


. -3,41 

— 


Note: e/k in °K; k 2 5 in — — - — ; AH in kcal/mole. 

1 For PE and PTFE k 25 is the solubility in amorphous phase; for 
FEP, k 25 is the experimental solubility. 

2 Reference 9. 

3 Reference 1. 

4 Heat of condensation. 

approximated by a linear relationship as has been observed also for other 
polymers. 16 A theoretical explanation has been given by Lawson. 17 

The activation energies and the diffusion constants at 25°C of PE , 
FEP, and PTFE are listed in Table 5. Rather surprisingly, the activation 
energies for the first two polymers are identical for all gases studied 
within the precision, of the measurements. According to the model by 
Brandt and Anysss\ s the activation , energy of -diffusion has two components; 
one -depends on the cohesive energy of the polymer and is proportional to 
the molecular diameter of the permeant t the other depends on the chain 
flexibility and -is proportional to the square of the diameter. Michaels 
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Table 5 


DIFFUSION DATA FOR DIVERSE GASES IN PE a FEP , AND PTFE 


Gas 

d 25 X 10 7 

e d 

PE 1 

FEP 

PTFE 1 

PE 1 

FEP 

PTFE 1 

n 2 

0.93 

0.95 

0.88 

9.0 

9.20 

7.12 

0 2 

1.70 

1.84 

1.52 

8.8 

8.29 

6.28 

co 2 

1.24 

1.05 

1.03 

8.5 

8.75 

6.84 

ch 4 

0 . 57 

0.30 

— 

10.4 

10.89 

-- 

c 2 h 6 

0.146 

0.047 

1 

12.5 

12.35 


c 3 h 8 

0.049 

0.0077 


13.6 

14.37 


c 2 h 4 

— 

0.098 

H 

— 

11.76 

— 


Note: D 2 5 in cm 2 sec; E^ in kcal/mole. 

1 Reference 9. 


and Bixler 9 estimate that for PE the contribution of the flexing energy 
to the activation ene.rgy is small, and accordingly find a linear relation- 
ship between the activation energies and the molecular diameters. In con- 
trast, the cohesive energy. of FEP is smaller than that of PE; but the FEP 
chain is stiffer. therefore, the energy associated with chain flexing 
should be significant in FEP, and the activation energies of diffusion 
should increase more than linearly with the molecular diameter of tne per- 
mesnts. 

Finally-, a comparison between the activation energies of diffusion 
in FEP and PTFE (Table 5) is of considerable interest. Although the co- 
hesive energy of PTFE is certainly quite similar to that of-FEP and Its 
chain stiffness is probably^higher , the apparent activation energi es for 
all the^three gases studies is about 2 kcal lower. Thus, easier diffusion 
pathways must be available -in the sintered PTFE ; they very likely are 
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micropores and/or grain boundaries. According to the manufacturer, commer- 
cial PTFE has a void volume of about 0,5^; nevertheless, the density of 
PTFE increases by several percent if compressed at 4 kbar/inch 2 9 Further 
study will be needed to show whether this increase is due to the removal 
of pores or to an increase in crystallinity. 

Conclusions 

Empirical correlations have been obtained in this and other studies 
between transport properties and molecular parameters of gases permeating 
through polymers. However ?™the apparently specific relationships observed 
for heat of solution and force constant, on the one hand, and activation 
energy of diffusion and molecular diameter on the other, may well be for- 
tuitous. For the gases studied, the molecular diameters are approximately 
proportional to the Lennard-Jones force constants (or their boiling tern - 
peratures) ; thus, one could have used either of these parameters and still 
obtained linear correlations for the heats of solution or of diffusion. 
Hence, a plot of versus E^ for FEP is approximately linear also (Fig, 
7); only C0 2 deviates again pronouncedly from the straight line. 

In the short discussion of the data we have pointed out the inability 
of the existing theories of gas-polymer interactions to explain satisfac- 
torily the experimental data for FEP and possibly for PE, aiso. It is 
beyond the scope of this present study to formulate models mo v e consistent 
with the observations. 

Finally, the conclusion that diffusion through PTFE proceeds at least 
in part along pores and gram boundaries may have a significance far beyond 
the aif fusion process. Other observed properties of bulk PTFE, which have 
been ascribed to the dense polymer structure itself, may also be due in 
part to the physical structure of the sintered material, A reevaluation 
of PTIE properties on the basis of this model may well be worthwhile. 
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LEGENDS TO FIGURES 


Fig. 1 Logarithm of the Permeation Coefficient P of Seven Gases in 
FEP, as a Function of the Reciprocal Temperature 1/T. 

Fig. 2 Logarithm of the Diffusion Coefficient D of Seven Gases in 
FEP, as a Function of the Reciprocal Temperature l/T. 

Fig. 3 Logarithm of the Solubility k 25 at 25°C of Seven Gases in 
FEP as a Function of the Lennard-Jones Force Constant e/k. 

Fig. 4 Heat of Solution AH g of Seven Gases in FEP, As a Function 
of the Lennard-Jones Force Constant e/k. 

Fig. 5 Logarithm of the Preexponential Factor k of Seven Gases in 

FEP as a Function of Heat of Solution AH°. 

S 

Fig. 6 Logarithm of the Preexponential Factor D of Seven Gases in 
FEP as a Function of Activation Energy o? Diffusion E^. 

Fig. 7 Heat of Solution AH as a Function of Activation Energy of 
Diffusion Ep. 
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Appendix C* 


VISCOUS FLOW AND DIFFUSION OF GASES THROUGH SMALL APERTURES 

by 

R. A. Pasternak and M. V. Christensen 
Polymer Chemistry Program 
Stanford Research Institute 
Menlo Park, California 94025 


ABSTRACT 

In recent studies of diffusion and permeation of gas*"; through poly- 
mer membranes, 1 ’ 2 the need was felt for a reliable method of determining 
detector sensitivities for a variety of gases. A simple calibration pro- 
cedure has been devised; it is based on. pressure controlled viscous flows 
of gases through .small apertures. This concept has been elaborated by a 
systematic study of the mechanism of gas transport through apertures; the 
results of this work are reported here. 


* To be published in J. Applied Physics, 
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Appendix C* 

VISCOUS FLOW AND DIFFUSION OF GASES THROUGH SMALL APERTURES 

by 

R. A. Pasternak and. M. V. Christensen 
Polymer Chemistry Program 
Stanford Research Institute 
Menlo Park, California 94025 

Experimental 

The experimental unit, shown schematically in Figure 1, is basically 
the polymer Permeation Analyzer,* which is used for measuring transport of 
gases and vapors through polymer membranes, 3 except that the membrane is 
replaced by a pinhole or a capillary. A carrier gas at constant flow rate 
sweeps past the aperture to a detector, in the present design a katharom- 
eter. The test gas, which is either the pure tracer gas or its mixture 
with carrier gas, flows past the upstream side of the aperture; gas mix- 
tures of known composition are prepared by adjusting the flow rates of the 
tracer and the carrier gases with needle valves (1 and 2) and measuring 
them individually with a soap bubble flow meter. 

A pressure difference, Ap , which . may be positive or negative . (i . e . , 
the pressure of the carrier gas is lower or higher than that of the test 
gas) is established across the aperture by throttling the two. gas streams 
with needle valves (3 and 4) located at the exit, ports; Ap is measured by 
a tube manometer, which is filled with a fluid of unit density. The pre- 
cision in the manometer reading is about 0.5 mm. (The pressure difference 
is here expressed in mm of the unit-density fluid; thus, one mm is equal 
to 98.1 dynes/'cm^ . 


* Manufactured by Dohrmann Instrument Co. (Inf otronics) , Mountain View ,* 
California. 
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In the present experiments, the carrier gas was helium, and the 
pressure in both gas streams was close to one atmosphere; the maximum 
pressure difference was less than 0,02 atmosphere. The temperature was 
maintained at about 25°C, 


Tracer gas passes through the aperture by two mechanisms, diffusion 
and volume flow. The rate of transport F of the tracer gas is related 
to the observed detector signal S T by 

F = k(f + V)S’ (1) 

o 

where k is the sensitivity of the detector for. the particular tracer gas, 

f is the flow rate of carrier gas* and V the volume flow of the test gas 
o 

through the aperture. At high pressure differences, when viscous flow is 
the principal transport mechanism, .the volume flow can be approximated by 
V ^ F / c 0 > where c^ is the partial pressure of the tracer in the test gas. 

Introduction of this expression into Eq. (1) leads to 

(f + V) = f 
o o 

This correction to f amounted at most to about 3<j 0 in the present study. 
The rate of transport can now be redefined by 

F = kS 

where S = (f Q + V)S T “ is the signal normalized for an effective carrier 
gas flow rate of unity,. 

Finally, it was established by variation in the flow rate f that 

_ o 

the detector, response was linear with concentration over the range 
studied. 



Theoretical 


Derivat ion of Equation for Total Flow of Tracer Gas Through Aperture 


Dilution and viscous flow occur simultaneously through the aperture, 
and their relative importance depends on pressure. One-dimensional total 
flow of tracer gas can be described by the differential equation 4 


dc 

dt 




(3a) 


where c is the partial pressure of tracer gas at any point X in the flow 
channel, t is the time, D the diffusion coefficient, and u the linear, 
convent ion al flow. rate of the test gas caused by external forces, here by 
a pressure gradient. 


At steady state, Eq. (3a) reduces to 


D 


a 2 c 

DX 2 


dc 

‘ 0 


The boundary conditions in our experiments are 


(3b) 


c = c at X = 0 
o 

c f* 0 at X = i 

Now, c q is the partial pressure of the tracer in the test gas, and i is 
the length of the channel. Integration (at constant pressure difference) 
between the limits given and rearrangement results in the expression for 
the rate of transport F of the tracer gas through the aperture 

(nr 2 ) uc 

o 

1-exp (-j^*) 


F 


(4a) 


where rrr 2 is the cross-sectional area of a cylindrical channel. In this 
integration, u and D are average values fo-r a particular configuration 
and are assumed to be independent of X and c. However, u is a function 
of the pressure difference A p; it can be positive or negative; in contrast, 

D should be independent of Ap> at least for the dilute gas mixtures, and 
should represent the interdiffusion coefficient of the tracer gas in the 
carrier gas. 

At A p = 0, the flow velocity u is zero, and Eq. (4a) reduces to the 
simple diffusion rate law 

rrr 2 Dc 

F = kS = - (4b) 

o o £ 

Introduction of F q and y = u^/D into Eq. (4a) leads to 

F_ _ S 1 ( . 

F q ~ S Q " l-exp(-y) 

For large positive values of y (say y > 4) , equivalent to high Ap values, 

Eq. (4a) reduces to 

F = kS = nr 2 uc (4d) 

o 

The parameter y is the ratio of the transport rates of the tracer 
gas by convectional flow and by diffusion; thus, in the intermediate 
pressure range the flux depends, on both -flow mechanisms. . The diffusion 
coefficient D can be determined from the signal at Ap = 0, and u can be 
evaluated as a function of Ap from the experimental data by means of 
Eq. (4c), and the value of D, 

B. Convectional Flow 

.In the preceding section it has been Shown how the -flow v.. Deity' u 
through the aperture can be_ derived from the .experimental, data. Since 
the total pressure in the experiments is high, approximately one atmosphere, 
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the flow can be expected to be viscous; the flow rate should be related 
to the pressure difference Ap across a cylindrical aperture by the modi- 
fied Poiseuille equation 


AP = pF U + 1.14 P u 2 


(5) 


where T| is the viscosity and p the density of the gas, and r and l are 
the radius and length of the aperture. 


The second term is an end correction which takes into account that 
the flow will be fully developed only after a finite distance. 

The equation holds if this term is less than about half the first term. 
Moreover, additional conditions must be satisfied : 5 


( 1 ) The fluid must be incompressible. Compressibility of the gas 
can be disregarded when the Mach number ,M (the ratio of flow and 
sound velocity) is smaller than one third. 

(2) The flow must be laminar. This condition holds when the Rey- 
nolds number R = rpu/Tl is smaller than 1,000. 

(3) The fluid .must be immobilized at the tube wall. 

Substitution of y = u^/D into Eq. ( 5 ) results in 


A 8 TI 0 D n 2 

Ap “ T 2 I ' y + 1# 14 p y 2 - + $y 2 (6a) 

For large values of y; i.e., at high pressure dif f erences, . the flow- vel- 
ocity can be expressed by Eq, (4d), ; and the experimental data should fit 
the equation 


Ap 
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= aS > bS 2 


( 6 b) 


Results and Discussion 


Two apertures were studied, a glass capillary and a pinhole in a 
metal foil* 


A . Capillar y 

The capillary radius r measured microscopically was 3.9 ±0.1 X 1CT 3 
cm, and its length measured with a micrometer, was 0.47 ± 0.005 cm. 

In Figures 2 and 3 the detector signal S is plotted versus the pres- 
sure difference £p across the capillary for .pure C0 2 , 0 2 , N 2 , CH 4 , C 2 H 6 , 
C 2 H 4 , and C 3 H g as tracer gases. Runs (not shown graphically) were also 
made with mixtures of tracer gases and helium. 

The data points for all series are described very well by. straight 
lines over the entire range of observation. It thus appears that the 
linear term of the Poiseuille equation (6). applies; i.e., that the 
correction term for end effects and the contribution of diffusion are 
both negligible. The highest flow velocity u obtained in any. of the 
experimental series was about 5 x 10 2 cm/sec; estimates show that at this 
velocity the two corrections are insignificant .and that the other condi- 
tions for the validity of the Poiseuille equation are satisfied also. 
Thus, the slope of the experimental line can be expressed in terms of the 
coefficients in Eq t (6b). 


dS 


c nr 
o 


dAp 8£T]k; 



(7) 


where A is a constant which can be. calculated from the dimensions of the 
aperture. 

Equation (7) „can be used -in two ways: If ' the viscosity of the fluid 

is known, the sensitivity of the detector .xan be derived; also the reverse 
if the sensitivity is known, the viscosity can be calculated. 

In Table I, column 2, the slopes obtained for the pure gases and for 
the gas mixtures studied are listed; they are expressed in mm/dynes cm“ 2 . 



Table I 

RELATIVE SENSITIVITY FACTORS OF KATHAROMETER FOR 
DIFFERENT GASES, AND VISCOSITIES OF 
GAS MIXTURES FLOWING THROUGH A CAPILLARY 


Gas 


dS/riAp 

(mm dyne" 1 cm 2 ) 

n 

(p, poise) 


o b s . 

■ 

C0 2 

100 

2.23 





100 

2.23 

148 

1.00 

1.00 


100 

2 . 25 




o 2 

100 

1.36 





100 

1 . 34 

202 

1.21 j 

1.17 

n 2 

100 

1.65 

175 

1.14 

1.09 

ch 4 

100 

2.28 

109 

1.32 

1.33 

C 8 H 6 

100 

3.75 





100 

3.77 

91 

0 .97 

0.94 

c 2 h 4 

100 

3.33 

101 

0.98 

0.98 

c 3 h 8 

100 

5.21 

80 

0.79 

0.72 

<M 

8 

17 .0 

1.68 

196 





8.6 

1.69 

196 

-- 

— 

0 2 

32.7 

1.26 

216 



—— 


20.4 

1.27 

214 

. — 

— 


13.1 

1.32 

206 


— 


9.1. 

1.34 

203 

— . 

— 


6.2 

1.34 

203 

— 

— 

c 2 h 6 

42.2 

3.26 

105 



_ _ 


21.0 

2.84 

120 

— 

. 


14.8. 

2.59 . 

132 

— 

— 


6.5 

2.20 

156 

— — 

— — 


Note: The viscosities of the pure gases are taken from 

the literature*. 6 those of the gas-helium mixtures 
are calculated from the experimental data.. 







Repeat runs are shown for three gases to indicate the reproducibility of 
the results. The viscosities for the pure gases are given in column 3; 
they are taken from the literature. 6 Column 4 shows the detector sensi- 
tivities k relative to that of C0 2 , and in column 5 the relative sensi- 
tivities derived from a tabulation by Kaiser 7 are listed. The agreement 
between the two sets of data is quite satisfactory. Moreover, the detector 
sensitivites for C0 2 , 0 2 , and N 2 were determined directly with gas mixtures 
of known composition. The relative values agreed with those given in col- 
umn 4 within the precision of the measurements. 

The absolute sensitivity of the detector can be calculated by means 
of Eq. (7) from the dimensions of the capillary. We obtain for C0 2 

f 

k^ = 4.82 x 10“ 7 cc/mm 

the direct calibration with the C0 2 '-He gas mixture gave 

k = 5.82 x HT 7 cc/mm 

This agreement is quite acceptable. The estimated uncertainty in the 
radius of the capillary was about 2.5$, and thus the experimental uncer- 
tainty in k^ is at least 10$, since the radius enters with its fourth 
power. Moreover, the Poiseuille equation is derived for an idealized 
case and the deviation of 20$ found here is easily possible for a real 
conductance. 

In the lower . half of Table 1 the data for gas mixtures are. given. 

Here, the viscosities are calculated from the slopes and the sensitivity 
factors which were derived for the pure gases. Figure. 4. is a graphical 
representation of the data. Obviously .the viscosities must approach at 
infinite dilution that of pure helium (196|j. poise); but we observe rather 
surprisingly, different concentration effects for the-three gases studied. 
However, dependence on . composition similar! to that for the He-0 2 mixtures 
has been reported for He-Ar mixtures,® A more extensive and .thorough 
experimental, and theoretical: study of the viscosity of gas mixtures would 
be worthwhile. 


B. 


Pinhole 


The pinhole dimensions were r := 5.0 ± 0.2 x 10 -4 cm and £ = 4 ± 1 x 
10- 4 cm. They were measured microscopically, but the length determination 
required cross-sectioning of the foil. 

The gases C0 2 , 0 2 , and N 2 (pure and diluted with He) were studied 
both for positive and negative pressure differences. 

The data were evaluated according to Eq. (4c); i.e., the parameter 
y = U£/D was computed from the ratio S/S q . A typical experimental run 
is shown in Figure 5 for a 22.2$ C0 2 ~He mixture. The parameter y is also 
plotted in the figure; within the precision of the measurements it is a 
linear function of Ap over the entire range of observation. Linear 
correlations of the same quality were obtained for all gas mixtures. For 
the pure gases, significant deviations from linearity were observed. This 
linear dependence of y on Ap indicates that the flow velocity u. obeys the 
one-term Poiseuille equation (5) with the viscosity independent of Ap, 
and that the diffusion coefficient D is constant also. 

The maximum flow velocity in these experiments was about 2 x 10 3 cm/sec. 
Both the Mach and the Reynolds numbers are at this velocity small enough to 
satisfy the conditions for viscous flow. However, for the mixtures the 
square term in Eq. (5) was calculated at the highest velocities to be 10 
to 20$ of the linear term. The experimental results. suggest that the 
second-order term is smaller than that estimated from the equation; or 
vnat it is compensated by other effects not considered here. Because, 
of the high density of the pure gases, however, the second term should 
be so large that the Poiseuille equation must break down completely;, still, 
approximately straight lines were obtained even for them, although the fit 
was less satisfactory. 

The viscosity and the. diffusion coefficient cannot be strictly 
constant over the whole pressure range studied, since the composition of 
the gas along the length of the aperture and its average must be a func- 
tion of-flow rate (or pressure difference). However, for the dilute 
gases, the concentration of the tracer being at most about 20$, the 


Oil 


r 


variation in concentration will not be large. Moreover, the diffusion 
coefficients of gas mixtures have been shown to bo barely dependent on 
composition, 9 and we have found the same for the viscosities of C0 2 -He 
and 0 2 -He mixtures flowing through a capillary (see above). For pure 
gases, however, both parameters must vary significantly with pressure 
since the composition changes from one extreme to the other. 

In Table II are shown the slopes of the y versus Ap plots, the 
diffusion coefficients derived from the signals at Ap - 0 by Eq. (4b), 
and the viscosities calculated from the slopes as defined by Eq. (6a) 


<SL_- £! . i- 

dAp 8 DTI 


The data for the pure gases are included also but are omitted from 
the discussion. 


Table II 


DIFFUSION AND VISCOSITY OF 
GASES FLOWING THROUGH A 10-p, PINHOLE 



Gas 

i 

dy/d Ap 

(cm 2 /dyne X IQ 4 ) 

D 

(cm 2 /sec) 

T1 

_(p, poise) 


C0 2 

100.0 

7.75 . 

0.123 

328 . 



22.2 

6.30 

0 . 189 

262 



16.0 

6.30 

0 . 189 

262 



11 . 2 

6.45 

0,180 

268 


0 2 

100.^0 

5.25 

0,158 

37G 



20.2 

4.82 

0.228 

254 



10.8 . 

4.98 

0.258 

243 


N 2 

100.0 

5.86 

0.163 

326 



~ 21.6 

5 . 25- 

0.228 

260 



18.3 

5.39 

7 0,226 

■ 256 




10.6 

5,37 

J) . 250 

232 _ 



9.1 

5.40 

. 0.252 

229 “ 
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The experimental uncertainties affect the parameters given in Table II 
to a varying extent. The slopes dy/dAp are obtained from the measured 
signal and pressure difference only and are therefore quite reliable. 

Table II shows that for a given tracer gas the slopes; and therefore the 
products T|D are approximately independent of the composition of the gas 
mixture. Thus; both Tj and D must be constant for a given tracer gas; 
as suggested already; the alternative explanation that they are inversely 
proportional to each other is excluded since they are both similarly re- 
lated to the cross sections of the gas molecules. 

The calculation of the diffusion coefficients involves the square of 
the radius and the length of the aperture; these parameters are associated 
with large experimental errors but are in common for all experimental ser- 
ies. However, the specific test gas composition, the measurement of which 
is not of highest precision,, enters individually into each run. Therefore, 
the variation of D with .composit ion is probably associated primarily with 
this experimental uncertainty. 

Finally the viscosity, by its mode of calculation, includes all the 
uncertainties associated with the diffusion coefficient. The pinhole ra- 
dius enters, however, with its fourth power into the calculations. 

In Table III the average observed diffusion coefficients and viscos- 
ities are compared with the interdiffusion coefficients calculated by the 
Stefan-Maxwell 1 0 formula, which usually agrees well with experimental data, 
a nd with the viscosity of helium taken from^the literature _____ - 

The observed values of D are smaller than the . calculated ones by ap- 
proximately a constant factor. This observed reduction in D is far too 
large to be explained by the experimental errors in the measurements of 
the pinhole dimensions but appears to be real; the diameter and length of 
the pinhole are of the same magnitude, and thus diffusional flow is im- 
peded by collision of the molecules with the walls. 

The viscosity of the dilute gases can be expected to approach that 
of pure He-, irrespective of the tracer gas. We find approximately_the 
same viscosity for the three dilute gases studied' The agreement with the 
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Table III 


COMPARISON OF OBSERVED, AVERAGED, AND 
CALCULATED DIFFUSION COEFFICIENTS AND 
VISCOSITIES OF DILUTE GASES IN HELIUM, PINHOLE DATA 


Gas 

D 

(cm 2 /sec) 

Dobs/D 

calc 

n 

( n- poise) 

ob.s 

calc 

obs 

He 

co 2 


0.44 

0.42 

264 


0 2 


0.62 

0.29 

264 

197 

n 2 


0.69 

0.41 

245 



literature value for He is surprisingly close ; the fourth power of the 
radius and the first power of the length of the aperture, both parameters 
being measured with low precision, enter into the calculation of the ex- 
perimental value. 
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DIFFUSION AND PERMEATION OF GASES THROUGH NITROSO RUBBER 

by 

R. A. Pasternak and M. V, Christensen 
Polymer Chemistry Program 
Stanford Research Institute 
Menlo Park, California 94025 


ABSTRACT 

Nitroso rubber, is of considerable interest for air space applications 
because of its resistance to strong oxidants and its good low- temperature 
flexibility. 1 A limited study of transport gases through this rubber has 
been carried out since such information is of importance for its potential 
use . 


* To be published in J. Applied Polymer Science. 
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Experimental 

Two nitroso rubber samples were studied. One sample membrane had 
been manufactured by Thiokol Chemical Corp. (CTVR-CTA cured) ; its thick- 
ness was 35.0 ± 0.5 mil. The second sample was prepared in our labora- 
tory from materials procured by Mr. D. Lawson from Jet Propulsion Labora- 
tory, Pasadena. Nitroso terpolymer (100 parts), Cab-O-Sil silica (20 
parts), and chromium trif luoroacetate (5 parts) were blended by standard 
procedures and cured at 307°F for 45 minutes; the resulting membrane had 
a thickness of 20 ± 3 mil. Both membranes were very soft and elastic. 

High purity (99$+) 0 2 , N 2 , C0 2 , CH 4 , C 2 H 6 , and C 3 H S from cylinders 
were used without drying. 

The instrument and the methods for measuring diffusion and permeation 
coefficients have been described elsewhere. 2 * 3 Determinations were made 
at about 60, 70, 80, and 90°C; for these thick. membranes , the. transmission 
rate* became; too small at lower temperatures to permit reliable measure- 
ments. 

In Table I the permeation and diffusion coefficients P and D at 60°C, 
and the heats of permeation and- diffusion E p and E^ are tabulated^ Since 
at most four experimental points were available for each Arrhenius plot 
(and only three for the thicker Thiokol membrane) , the uncertainties in 
the heats, which are proportional to the slopes, are quite large, possibly 
0.5 kcal/mole . The ^solubilities k, which are included in the table, are 
derived from the permeation and diffusion data in the usual way , k = P /D . 


Tabic 1 


TRANSPORT DATA FOR NITROSO RUBBER 


(a) SRI Membrane 














The data obtained for the two membranes agree quite closely, if one takes 
into account the rather large uncertainties in the measurements. 


The transport parameters follow a Systematic pattern. The diffusion 
coefficients decrease and the solubilities increase with increasing boil- 
ing point (or molecular diameter) of the gases; the inverse holds for the 
heats of diffusion and solution. Similar correlations have been found for 
polyethylene 4 and FEP teflon . 5 

The transport properties of N 2 , 0 2 , C0 2 , and CH 4 in nitroso rubber 
(extrapolated to 25°C) and in natural rubber , 6 ’ 7 are compared in Table II 
The diffusion coefficients and the activation energies of diffusion of the 
gases (except of CO ) are about the same in the two rubbers, but the solu- 
bilities are 2 to 4 times larger in the nitroso rubber, which is therefore 
more permeable. The somewhat different pattern for C0 2 may indicate strong 
specific interactions between C0 2 and the nitroso rubber. 


* 


Table II 

COMPARISON BETWEEN NITROSO (Ni) AND NATURAL (Na) RUBBER 


Gas 

d 25 X 10 6 

- E r 

) 

' K 2 5 


.. Ni 

Na 

Ni 

Na 

Ni 

Na 

Ni 

-Na - 

N* 

0.92 

■ 

6.9 

8.7 

0 J23 

0.005 

2.1 

0.6 

0 2 

1.35 

■ 

6 .3 

8.3 

0.46 

0 . 112 - 

6.2 

1.8 

co 2 

0.25 

U 

8.8 

8.9 

5.8 

0.90 

in 

H 

9.9 

ch 4 

0.86 

0.1 

7.5 

8.7 

0.43 

0.25 

2.8 

CM 

CM 


Note: 


D in cm* /sec, 


k 


cc (STP) 
cc atm * 


r, CC (STP) cm 

P m — 5 — . 

cnr sec atm 


E^ in kcal/mole. 
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Through Plastics by R. A. Pasternak and James A. McNulty, 

Modern Packaging Magazine, April 1970. 
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